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DC CALCULUS 


LUIGI AMBROSIO AND JEROME BERTRAND 


Abstract. In this paper, we extend the DC Calculus introduced by Perelman on finite dimen¬ 
sional Alexandrov spaces with curvature bounded below. Among other things, our results allow 
us to define the Hessian and the Laplacian of DC functions (including distance functions as a 
particular instance) as a measure-valued tensor and a Radon measure respectively. We show that 
these objects share various properties with their analogues on smooth Riemannian manifolds. 
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1. Introduction 

In this paper, we investigate the DC Calculus introduced by Perelman in HZ] on the manifold 
part X* of a finite dimensional Alexandrov space X with curvature bounded below. The term ^^DCP' 
stands for difference of concave functions. The main point is that, contrary to the set of (semi- 
)concave functions, the set of DC functions on Euclidean space is an algebra. Perelman showed 
the existence of an atlas on A* compatible with the distance whose transition maps are Euclidean 
DC functions. Moreover, he also showed that DC functions on Alexandrov spaces (defined in 
geometric terms involving geodesics) are standard Euclidean DC functions when read in a chart. 
This notion is also natural if one considers the particular case of closed convex hypersurfaces in 
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Euclidean space where the standard charts are inverse of convex functions. The DC structure is 
also considered in na in connection with properties of the heat kernel on Alexandrov space. 

Besides providing a complete proof of the results related to this calculus discussed in m , we 
also generalize them in some aspects. Note that we choose to present our results in the case of ten¬ 
sors (mainly covariant ones) instead of differential forms as it is in the case in the aforementioned 
papers. However, it is then a standard algebraic matter to infer from our results the corresponding 
statements on differential forms (up to considering the orientable double cover of A* if necessary). 
More precisely, the calculus we develop allows us to consider general bounded BV functions (com¬ 
pared to standard BV functions, the BV functions satisfy an additional weak continuity property, 
see ®), thus removing the restriction to bounded BV functions with no jump part in their de¬ 
rivative (these functions are called BVq functions in Perelman’s paper )0 as in nznn]. The main 
result on tensors can be summarized as the existence of a covariant derivative operator on tensors 
S with BV components (thus DS is a measure-valued tensor) that satisfies the same properties 
as on a smooth Riemannian manifold. We are also able to define the norm of a measure-valued 
tensor that coincides with the standard norm in the case of function-valued tensor and gives us, 
for example, an intrinsic notion of total variation \df\ for a function of bounded variation / on X*. 

Our particular interest in covariant tensors stems from our will to define the Hessian of a distance 
function dp as a measure-valued tensor. Indeed, even in the case of a smooth Riemannian manifold, 
the distributional Hessian Hess dp of dp is not, in general, absolutely continuous with respect to 
the volume measure, see for instance m for a nice and comprehensive discussion of this fact. The 
extra term in Hess dp is of jump type, namely it is concentrated on the Cut Locus of p which is 
(TV — l)-dimensioiial in general, being N the dimension of the Riemannian manifold. For example, 
this phenomenon arises on the real projective space, details can be found in m- Consequently, 
even on a smooth Riemannian manifold, the Riemannian gradient Vdp does not belong to the class 
BVq defined by Perehnan in general. Our results allow us to define the Hessian of any distance 
function dp as a particular instance of a result for arbitrary DC functions. This generalization 
requires fine properties of BV functions and a thorough study of the jump part of the Hessian. 
Technically speaking, the primary issue is that, given two BV functions /, h defined on an open 
set of Euclidean space, the standard chain rule 

d{fh) ^ dh ^ df 

dxi dxi dxi 

does not hold true when /, h have non-trivial jump parts in their derivatives (see however Lemma [4.8ll 
However, even when jump parts in the derivatives are present, it may happen that cancellation 
properties occur, so that some standard formulas and definitions of tensor calculus retain their 
validity (see Proposition 15.31 and Proposition 16.21 as a particular instances of this phenomenon). 

As a consequence of this extended DC Calculus, we also prove that the DC Laplacian (namely 
the trace of Hessian with respect to the Riemannian metric) of DC functions coincides with the 
weak Laplacian defined through integration by parts (see Proposition l5.9|l . This fact has been used 
in |2fl| without an explicit proof. 

We are also able to prove that for a DC function / and BV vector fields X, Y, the Hessian of 
/ is related to the covariant derivative in the same way as in the smooth case in the sense that 

Hess/(A,y) = D{df{Y)){X) - df{DxY) 

(where D stands either for the differential of function or the (measure-valued) covariant derivative 
of vector field, depending on the context). We also prove that the classical formula 

= iv(|ViA|2) 

remains valid in this setting. 


^We invite the reader to compare our Section [31 with lemma p5 in Perelman’s paper 1171 . and our Section 14.41 
with Section 4.3 in the same paper 
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The formula above is then used to prove that the Hessian we define satisfies an integration by 
parts formula in the same spirit as in the r 2 calculus [5], see Proposition 16.21 This consistency 
result is particularly relevant in connection with the coordinate-free approach to calculus in metric 
measure spaces, which works particularly well under Ricci lower bounds on curvature, still expressed 
in terms of the r 2 tensor, see m- We know after (SUI and [Hj that spaces with sectional curvature 
bounded from below satisfy Ricci lower bounds in this sense, see also [55] for a very nice application 
of these calculus tools to Alexandrov spaces. Nevertheless the connection between this viewpoint 
and the one developed in this paper is not yet completely clear, since the calculus developed in 
[El El allows to handle, so to speak, only sections of the tangent bundle defined up to ^negligible 
sets, while the calculus in coordinates allows to handle sections of the tangent bundle defined up 
to {N — l)-dimensional sets. 

The DC Calculus allows Perelman to improve an important result of Otsu and Shioya in m 
where the authors prove that the distance on Alexandrov space X derives from a continuous (in a 
weak sense, see Theorem B in |16| i Riemannian metric defined almost everywhere on A*. Perelman 
proves that the Riemannian metric components, when read in a chart, are actually functions of 
local bounded variation. In |3], we improve these results in the case of surfaces. 

This paper is organized as follows. In the first part, we introduce the function spaces related 
to the DC Calculus. We then introduce the notion of DCq Riemannian manifold, which provides 
a natural setting where the DC Calculus can be defined. In the subsequent part, we study the 
covariant tensors and their covariant derivatives on a DCq Riemannian manifold. Due to the low 
regularity, we use the old-fashioned approach consisting in defining tensors in charts and imposing 
a compatibility condition (see |21| for more on this point of view). The antepenultimate section is 
devoted to the Laplacian and Hessian of a DC function. The penultimate one is devoted to the 
integration by parts formula involving Hessian. In the last part, we establish that the subset X* 
of a finite dimensional Alexandrov space X is indeed a DCq Riemannian manifold. 

2. Function spaces for DC Calculus 

2.1. DC, BV and Cw functions. 

Definition 2.1 (DC functions). Let D C be an open set. A function / : D —>■ R is said to be 
a DC function on D if it is locally representable in D as a difference of two semiconcave functions. 

A similar definition will be used for vector-valued maps, arguing componentwise. Since SC 
functions are defined by a local property, a simple partition of unity shows that DC functions in 
an open set can also be globally represented by the difference of two SC functions. It is also easy 
to check that the class of DC functions is a vector space. In the note na, Hartman proved that 
the space of DC functions is indeed an algebra, i.e. stable under multiplication. A more general 
statement, that encompasses also this property, is given below. 

Theorem 2.2 (Stability of DC). Given f : U C —> V C R'^ and g : V —> R with U and 

V two open sets, the function g o f is DC on U if f and g are DC on U and V respectively. As a 
consequence, the set of real-valued DC functions on U is an algebra. 

It is a classical result that the second distributional derivative of a SC function in an open 
set V C R^ is a symmetric matrix-valued Radon measure in V whose positive part is absolutely 
continuous with respect to the Lebesgue measure , and that SC functions admit at almost 
every point a second order Taylor expansion [T]. In particular, by linearity, the first derivatives of 
a DC function are functions of locally bounded variation. For later use, let us quickly review some 
basic definitions and results on BV functions on Euclidean spaces. 

Definition 2.3 {BV functions). Let V be an open set of R^. A function / : V — > R belongs 
to BV{V) if / G L^{V) and the derivative in the sense of distributions of / is representable by a 
vector-valued measure 


D/= (Di/,...,D^/) 
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with finite total variation in V. Equivalently, \Df\{V) < oo and 



Vi = 1,... ,iV, (/) G Lipc{V). 


Here Lipc{V) stands for the set of Lipschitz functions with compact support on V. 

The function / is said to be in BViodV) if it is a BV function on every open set V (s V. 


Again, a similar definition can be given componentwise for maps / : V —> and in this 

case Df will be viewed as a (M x Af)-matrix of measures Denoting by k > 0, the 

fc-dimensional Hausdorff measure (the ambient space is irrelevant and it will not appear in the 
notation), we recall that a set E C R.^ is said to be a-finite w.r.t. if E is the union of 

an increasing sequence of sets with finite measure. If /c > 1 is an integer, we say that a 
set E C R.^ is countably -rectifiable if there exist countably many closed sets Ci C R* and 
Lipschitz functions fi : Ci —>■ R^ such that 

OO 

jt^dE\\jMCd) = o. 

i=0 


Remark 2.4. Throughout this paper, (•,•) stands for the standard Euclidean inner product and 
II • II 2 for its induced norm. 


Given a locally integrable function / we say that x is an approximate continuity point of / if 
there exists f{x) G R satisfying 

(1) lim'f \f{y)- f{x)\dJ^'^{y) =0, 

where, here and in the sequel, Bp(x) denotes the open ball with radius p and center x and jg f dp 
denotes the averaged integral, i.e. Jg f dp/p{B). The complement of the set of approximate 
continuity points will be denoted hy Sf. We shall also use the same notation and concept for 
vector-valued functions, arguing componentwise. 

Theorem 2.5 (Decomposition of Df). Let V C R^ be an open set. For all f G BViofiV) the 
following properties hold: 

(a) |D/| vanishes on -negligible sets. 

(b) Df admits the mutually singular deeomposition 

(2) Df = D°-^f + D^d = D^‘'f + 

where D‘^‘^f is the absolutely continuous part w.r.t. Lebesgue measure , D^'^f is the 
singular part, whieh can be further split in jump part D^^ f and Cantor part D'^°' f, the 
former coneentrated on a set a-finite w.r.t. and the latter vanishing on all sets 

with finite measure. 

(c) D^^ f is eoncentrated on the set of approximate jump points Jf and the set Jf is a eount- 

ably J^^~^-reetifiable set. Moreover up to a -negligible set, Jf coincides with the 

approximate discontinuity set Sf of f. 

(d) By the very definition of Jf, all points x G Jf one can identify two values f^{x), the 

so-called one-sided approximate limits, such that an approximate continuity property holds 
by averaging along an half space passing through x. More preeisely, there exists a Borel 
function v{x) G such that 

lim/ \f{y)-f^{x)\d^^iy) = 0, 

where 

H+(x,u(x)) = {y G Bp{x)-, {y-x,v{x)) > 0}, B/;{x,v{x)) = {y G Bp{x);{y - x,v{x)) < 0}. 
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The triple (u(x),/^(a:), / (x)) is uniquely determined up to v{x) i—>■ —v{x) that maps 
f^ix) to f~{x) and vice versa. The jump part of Df is equal to 

- /-(x))u(x)^^-iL J/. 

For a proof of this proposition, we refer to the the book [U Chapter 3] (see also Dehnition 12.81 
below). 

We will also use the following result. 

Theorem 2.6. Let F : U ^ V be a biLipschitz homeomorphism between open subsets ofM.^ such 
that, either det dF’ > 0 -a.e. on U or det dF < 0 Jff^-a.e. on U. Then, for all f G BV(U) 

the function f o F~^ belongs to BV{V) and 

\D{foF-^)\<iL^pF)^-^F^{\Df\), 

where stands for the pushforward of a measure (see Section \3.S\) . 

For a proof, we refer to [H Theorem 3.16]. 

Lemma 2.7. Let f G BV{Ll), with C open and connected. If D^^f = 0 and, for any ball 
Brix) (E 11 one has |/|~^ G L°°{Br{x)), then the sign of f is constant. 

Proof. By the connectivity of fl, it suffices to show that if i? is a ball, / G BV{B), |/|“^ G L°°{B) 
and = 0 in B imply that / has constant sign in B. Notice that D^'^f = 0 in B implies 

n Sf) =0. If / has not a constant sign, then both {/ > 0} and {/ < 0} have positive 
.^-^-measure. By the coarea formula we can hnd s, t > 0 with |/| > max{s,l} in B such that 
{/ < —s} and {/ > t} have finite perimeter. By our choice of s and t and by the hypothesis on 
1/1“^, B is partitioned up to -negiigihle sets in two sets of finite perimeter which therefore 
have in B the same essential boundary (namely the set of points of density neither 0 nor 1). On 
the other hand, since o implies that all sets 

{yGBp{x): \f{y)-J{x)\>e] e>0 

have 0 density at x whenever x ^ Sf, it follows that any point in the intersection of the essential 
boundaries of {/ < —s} and {/ > t} belongs to Sf. It follows that the essential boundary of these 
sets is ^-negligible in B, and this forces one of the two sets (by the relative isoperimetric 

inequality for sets of hnite perimeter) to be J^^-negligible and the other one to coincide, up to 
J^^-negligible sets, with B. □ 

Definition 2.8 (Precise representative of a BV function). Let H C open and / G BV(fl). 
Recall that, by definition, at all points x G fl \ Sf there exists a number, called approximate limit 
and denoted f(x), satisfying /g \f — f{x)\ dy = o{r^). In addition, Lebesgue’s theorem gives 
that f = f -Sf'^-a.e. in LI. 

Analogously, at all points x G Jf one can identify two values f^{x), the so-called one-sided 
approximate limits, where an analogous property holds by averaging along an halfspace passing 
through X (see Theorem 12.51 ('dll. The precise representative /* is defined on Ll \ {Sf \ Jf) by 

ifxGn\5'/; 

+ f~{x)) iixGJf. 

and it is left undefined on the .;?f'^“^-negligible set Sf \ Jf (in particular, Df = Df*). In the rest 
of this paper a BV function / will be always identified with its precise representative. 

We refer to [d] Corollary 3.80] for more on the precise representative. Notice that the use of the 
representative is necessary in several results in this paper. For instance, the precise representative 
f* will be used when the chain rule formula [U Theorem 3.96] is applied in our setting, see also 
Lemma 14.81 
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2.2. Cw functions, geometric BV functions and Radon measures. 

Definition 2.9 (C^, and Cw,o functions). Let C K.^ be an open set and / : —>■ K. We say 

that / G if there exists a set S, cr-finite with respect to such that: 

(a) f\n\s is continuous; 

(b) /|a\S is locally bounded in fl, i.e. for all a; G one has |/| < oo for some ball 

Br(x) (E f2. 

We say that / G Cu,^o{^) if, in addition, S can be chosen to be -negligible. 

Notice that (a) does not imply (b), unless S = %. Occasionally we will use the notation (7^(0, S) 
or Cw,o{^, S) to emphasize the role of the set S in the defining properties of the classes (7^,(0) 
and Cw,o{^)- 

It is obvious that if / G S), then Sf C S and / = / = /* on O \ 5. In particular, for 

/ G BVioci^) 7 tbe precise representative coincides with / 7^^“^-a.e. on O. Most of the 

results stated below in the class C^i^) O BV{il) (for instance Lemma [4.81 belowl could actually 
be stated and proved in the class BV (O) using the precise representative; however the use of the 
spaces (7m (O) and (7m,o(f^) is compatible with our geometric applications and it simplifies some 
technical aspects. 

Notice also that the Euclidean gradient V/ of a DC function / : fl C —>■ K. (arbitrarily 

defined on the set of points where / is not differentiable) belongs to (7^(12). The following definition, 
then, is natural. 

Definition 2.10 (The class DCo{fl)). Let 12 C be an open set. DCo{Q) is the collection of 
all functions / G DC(fl) such that V/ G ((7m,o(f2))^. Analogously, Z2(7o(I2, 5”) is the collection of 
all functions / G DCiyi) such that V/ G ((7m_o(f2, 5*))^. 

For later use, let us also introduce 

(3) SV(f2) := BViociS^) n Cm(I2) and SVo(I2) := BVioci^) n (7m.o(I2). 

Analogous definitions can be given for the subclasses SV(I2, S), SVo(I2, S). 

The natural Radon measures related to these spaces are introduced in the following definition. 

Definition 2.11 {QM.{Vl) and QM.q{VI) Radon measures). The set of Radon measures on 12 that 
vanish on 7f^^“^-negligible sets is denoted by QM.{D), whereas the set of Radon measures on 12 
that vanish on 7^^“^-finite sets is denoted by t/Alo(f2). 

Remark 2.12. In our terminology, a Radon measure p on 12 has finite total variation only on 
compact subset of 12. Consequently, if p is a signed Radon measure, p(i?) makes sense only for 
Borel sets B with compact support in 12. Note that the derivative of / G SVo(12) has no jump 
part, therefore according to Theorem 12.51 it belongs to QA4o{^) . 

3. DCo Riemannian manifolds 

3.1. Definition and properties of DCq Riemannian manifolds. We start with the definition 
of the class of manifolds we are interested in. 

Definition 3.1 {DCq Riemannian manifolds). A countable at infinity N dimensional topological 
manifold X* is said to be a DCq (open) Riemannian manifold if the following two properties hold: 
(a) There exist a singular set S C X* and a maximal atlas of X* {Ua, (/'a)aGA made of charts 
cj)a : Ua ^ which are homeomorphisms onto their image such that for any a, (3 G A, 

4>aiUa n S) is 7^^“^-negligible, 

and any transition map 

F = (j)pO : Ua ■■= (jiaiUa H Ujj) ->• (j)p{Ua C Up) =: Up 
belongs to DCQ{4>a(Ua H Up),4>a{Ua n S)). 
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(b) X* is endowed with a Riemaiinian metric g defined on X* \ S whose components gij, when 
read in a chart {U, cj)), belong to BVo{4>{U),4>{U D 5)) and satisfy 

(4) -r-Tlb|| 2 > 9 ij{x)p^Pj > c{x)\\p\\l for all p e for all X e (/i(C7\ 5), 

with c and c~^ locally bounded in 

Remark 3.2 (Sign of determinant of transition maps). The determinant det dF of the transition 
map is a BV function with no jump part of derivative. Therefore, Lemma 12.71 and condition 
(b) in Definition 13.11 give that the sign of det dF is '^-almost everywhere) constant on the 
connected components of the domain of F. This allows us the application of Theorem l2.61 namely 
the invariance of the BV property under composition with F~^. 

Remark 3.3. The Riemannian metric in (b) is in particular a covariant 2-tensor. This means that 
the components of g when read in two different charts satisfy a compatibility condition. This 
condition is recalled in Section 01 

Remark 3.4. Note that a DCq Riemannian manifold is in particular a Lipschitz manifold. There¬ 
fore, there is a well-defined first order calculus for locally Lipschitz functions. Namely, a function / 
is said to be locally Lipschitz ii fo(j)~^ is a locally Lipschitz function with respect to the Euclidean 
norm when read in an arbitrary chart (U, (/>). Throughout this paper, “locally Lipschitz” will always 
refer to this definition. 

We can extend all the definitions of Section 12.21 to functions or Radon measures on a DCq 
Riemannian manifold, by requiring that the corresponding definition given there is satisfied when 
the object is read in any chart ([/, (j>). Namely, a function / belongs to the space “X” if and only 
if f o(j)~^ belongs to the space “X” in the Euclidean sense, and a Radon measure p belongs to the 
space “X” if (j)^p belongs to the same space in the Euclidean sense. The well-posedness of these 
definitions is guaranteed by the fact that a DCq transformation map is in particular a biLipschitz 
homeomorphism, and the fact that sets which are cr-finite or negligible w.r.t. are invariant 

under biLipschitz maps. 

3.2. Properties of the metric. We end this part with technical results that will be used in the 
sequel. 

Lemma 3.5. Let X* be a DCq Riemannian manifold. We set G~^{x) = {g’’^{x))i<ij<i\[ the 
inverse at x G (f>{U\S) of the matrix {gij{x))i<ij<N read in a chart {U,(j)). Then, G~^ belongs to 

svo(</.(c/))^b 

Proof. Let D = ip{U) and let S be the -negligible singular set oi g, in these local coordinates. 

Using the expression of the inverse of a matrix in terms of the matrix of cofactors and the fact 
that BVq is an algebra, it suffices to prove that l/'g, with 'g = det{gij), belongs to BViodLt) and is 
locally bounded in D, according to Definition 12.91 Local boundedness of l/'g follows immediately 
by condition (b) in Definition 13.11 The fact that l/'g belongs to BVioc{Ll) follows from the local 
boundedness and from the chain rule formula for the left composition with a Lipschitz function 
(see for instance [H Theorem 3.96]), which gives 

D— = - nDg 

9 t 

provided we work with the precise representative (recall also that g has no jump part in its derivative 
since it belongs to BVo{Ll)). □ 

For later use, we also mention the following result whose proof is along the same lines as the 
one above, thus is left to the reader. 
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Lemma 3.6. Let X* he a DCq Riemannian manifold and let G be the BVo-metric read in a chart. 
Then, understanding the derivatives in the sense of distributions, one has 



To conclude this part, we notice that as in the smooth case, a Riemannian metric can be defined 
locally first and then be extended to a global tensor on the DCq manifold. Indeed, it is easy to 
check that Q is preserved by DCq transformation map and the lemma below guarantees the 
existence of DCq partition of unity. 

Lemma 3.7. Given a locally finite atlas (Ui, on X*, there exists a partition of unity 

subordinate to where ipi are compactly supported Lipschitz DGq functions in Ui. 

Proof. X* is paracompact, thus there exists a locally finite subcover {U[)i^i of {Ui)i^i with (s 
Ui. Since in Euclidean spaces one can always find, given open sets A (s B, a. function (f G Cf°{B) 
with 0 < (/) < 1 and ^ = 1 in a neighbourhood of A, we can pull these maps (with A = Ai = (fifU^), 
B = Bi = (j)i{Ui)) to obtain DGq functions ipi and build out of them the desired partition of 
unity. □ 

3.3. Pull-back of function and measure on a DGq Riemannian manifold. Recall that, for 
a proper map F : X — > Y, the push-forward operator maps Radon measures in X to Radon 
measures in Y via the formula F^fi{B) = p,{F~^{B)) for any Borel set B with compact support in 
Y. Equivalently, since Borel functions can be approximated with simple Borel functions, one can 
characterize Ej/i via the change of variables formula 



In particular, we obtain the useful formula 

(6) Fiiik o F)fi) = kF^n for any locally bounded Borel function k : Y —>■ R, 

since (by applying the chance of variables formula with y? = kxs) 



for any Borel set B with compact support. 

For geometrical purposes, we need also to introduce the pull-back of measures, in a form that 
takes into account the Jacobian determinant of F, see © below. 


Definition 3.8 (Pull-back of functions and measures). Let X* be a DGq manifold. Let F : U 
R be a DGq transition map. Given / G Gw{V) we define the pull-back of / through F by 


F*{f) = foF. 

Analogously, given p, G QXi{V), we define the pull-back of p through F by 


{F*{p),'ijj) = if o F ^IdetdE ^\p{dx), 



where if is any compactly supported and bounded Borel function on U. Equivalently, taking the 
change of variables formula for F^ into account, one can write 


(8) 


( 7 ) 


I F*{p) = {F-^)i{\detdF-^\p). 

Combining ([0]) and © we immediately get a formula “dual” to (O , namely 
I F*(kp) = {k o F)F*{p) for any locally bounded Borel function k. 
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If ^ and '0 is any compactly supported and bounded Borel function, the change of 

variable formula of for blLipschitz homeomorphisms gives 

J ijj{x)dF*{p,){x) = j 'ijj{F~^{x))\deidF~^\p{x)dx 
= J ip{x) p{F{x)) dx. 

Thus, 

(9) F*{pjF^) = {F*p)Jf^ 

and the two definitions given above, for functions and measures, are mutually consistent. 

Remark 3.9. Note that x i—>■ Idetda,!^”^! belongs to BVo(V) thus, in particular, it is defined 
/r-almost everywhere whenever p £ G.M(V). 

3.3.1. Approximation result and geometrical consequences. In this part we prove a proposition 
whose corollary will be used many times in the rest of the paper. 

Proposition 3.10. Let 12 be an open subset o/M^ and h : 17 —>■ K. bounded. Assume that either 
h £ Cw,o(Gl) and p £ t/AI(17), or h £ (7^(17) and p £ GMo(Ll). Let (/ie)e>o be Radon measures 
with finite total variation in 17 absolutely continuous with respect to such that p^ ^ p in the 
duality with C'c(17). We further assume that 


( 10 ) 


limsup |p,£|(17) < |/r|(17) < oo. 

e4-0 


Then hp^ —>■ hp in the duality with C'{,(17), the class of bounded continuous functions in 17. 


Proof. The leading idea of the proof is that the set of discontinuity points of a suitable restriction 
of h is negligible with respect to p and p^. From [U Theorem 2.2], we infer that \pg\ —^ |^| in the 
duality with Cb{Ll). As a consequence, we clearly have that pf —>■ p^ in the duality with Cc)!!). 
Combining this remark with the analogous decomposition of h into its positive and negative part, 
the boundedness of h reduces the proof to the case when Pe > 0 and h is nonnegative and bounded. 

Let us consider the case h £ Cw{Lt, S) and p £ SVo(17), the proof in the other case is similar. 
We then define the lower semi-continuous envelope /i_ : 17 —>■ M 

h-{x) = infjliminf h{xn)', Xn £ ^ \ S and lim Xn = x}. 


We define the upper-semicontinuous h.^. envelope analogously. Now, using that /i_ > 0 is lower- 
semicontinuous and p > 0, we get (see for instance [H Proposition 1.62] for a proof) for any 
nonnegative continuous function if 


lim inf 
e4,0 


J Tph- dpe > J 


0/i_ dp. 


Analogously, if ip is also bounded, we have 




lim sup / iph+ dpe < 

eX-Q 


iph+ dp. 


Now, the continuity assumption on h yields /i-|_ = on 17 \ S, so that 


0/l_|_ dpe = / iph- dps = / iph dpe 


and 


J 0h_|_ dp = J Iph- dp = J Iph dp. 


This proves the convergence of hp^ to hp in the duality with bounded and nonnegative continuous 
functions. The general case can be achieved splitting the test function ip in positive and negative 
part. □ 
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Corollary 3.11. Let F : U —> V be a DCq transition map. Assuming that h G BV(V), the 
following chain rule in the sense of Radon measures holds: 



for any i G {!,■ ■ ■ , N} and F = {Fi, ■ ■ ■ , Fn). 


Proof. We set = h* with p^ a family of mollifiers. Then, hg and h^o F are locally Lipschitz 
functions. Thus, (HU holds with he instead of h. Since hg o F ^ h o F in one has 

a(k OF) , = 

oxi axi 

in the duality with C^{U), and then in the duality with Cc{U). On the other hand, using 
F*(hp) = {ho F) F*{p), we can rewrite the right-hand side of (fTTl) as 



Now, pe = dhe/dysM’^ satisfy the hypotheses in Proposition [3T0] (a proof is given in [H Theorem 
2.2] for instance), dh/dys G QM.{V), \ (ietdF~^\dFs/ dxi o F~^ G o{V) therefore Proposition 
ISini yields 


n I? 


dh 

dys 


in the duality with Cc{U). This implies, using F*{p) = F^ ^(| detcif ^\p), 

F’{^oF-^fa.^F}^F-i^oF-‘ 

UXi ^ys (^Xi ^ys 


in the duality with Cc{U). By combining this convergence with ([ 5 ]), we get (ED- 


□ 


3.3.2. Measure induced by a system of measures. 

Definition 3.12 (System of Radon measures). Let X* be a DCq Riemannian manifold with charts 
{Ua, 4>a)a&A- A family of Radon measures {pa)aeA, where pa is a Radon measure on 4)a{Ua) is 
called a system of Radon measures if it further satisfies the compatibility condition F*{pp) = pa 
in (faiUa n Up). 

Note that the above condition is not equivalent to F^{pa) = pp. However, we have the following 
result. 

Lemma 3.13. Let F : U —> V be a DCq transition map. Suppose we are given p G GA4{V) and 
set p = F*{p) G QM.{U). Then, 

F^{V det Gp) = V det G p, 

where G and G are the metrics in the respective coordinate systems. 

Proof. Let if he a bounded Borel function with compact support. In this proof, we will need the 
fact that 

G = *{dF){GoF) {dF) 

and 

G-i = (d^(.)R-i) (G-ioF)‘(d^(.)F-i). 

where stands for matrix transposition. These relations are nothing but a reformulation in terms 
of matrices of the compatibility condition applied to the covariant 2-tensor g in different charts. 
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This compatibility condition is recalled in Definition 14.41 We infer from the equalities above that 
\/det G = Vdet G o F~^ \ det dF~^\. Thus, we get 


(Fjj ^{\/detGn),tp) = j ipoF \y)\/det G(y)/r((iy) 

= J Ip o F~^{y)\/det G o F~^{y)\ det dF~^\{y)'jl{dy) 
= (Vd^F*(M),V>), 


where we used the identity F*(hy,) = {ho F)F*{y) to get the last equality. 


□ 


For later use, let us also point out that combining the fact that F is a homeomorphism, F*(/i) = 
detdF“^|^), and | detdF“^| > 0 yields 

(12) \h\=F*{\y\) 

by uniqueness of the Radon-Nikodym decomposition (see [H Corollary 1.29] for a precise state¬ 
ment). 

Lemma 3.14. Let X* be a DCq Riemannian manifold with charts {Ua, 4>a)aeA- Then, any system 
of Radon measures (/ia)aeA with fj,a S QM.{4>a{Ua)) for all a G A, induces a Radon measure 
fj. G GA4{X*), characterized by 

(13) = / \/det Ga (x) dyLa{x) for any Borel set B C Ua- 

Urb) 

The total variation of this measure is then given by 

(14) \y\g{B) = / v'd^ a{x) d\yLa\{x) for any Borel set B C Ua- 

Urb) 

Morever, if fia G GMoi(pa{Ua)) for all a G A, then y, G 0Aio(Ai*). 


Proof The proof is based on a standard argument involving partitions of unity. Recall that, 
according to Lemma 13.71 there exists a locally finite partition of unity {'ipa)a&A subordinate to 
{Ua)aeA where ipa are compactly supported Lipschitz DCq functions. Then, given any bounded 
Borel function y with compact support, we set 

/ Xdy-'^ i^aX) ° ^ {x) \J det (x) dy^ {x) ■ 

d X* a^A^’t>a{Ua) 

Using Lemma [3.131 and (TT^ . it is a standard fact that the measure is well-defined (i.e. that it 
does not depend on the partition of unity) and that it satisfies (IT^ and (11411 . We refer to El 
or to m for more details. Moreover, it is immediately seen that y G GA4o(X*) whenever all 
yaGQMo((pa{Ua))- □ 


4. Tensors on DCq Riemannian manifolds 

Our goal in this section is to prove that there is a well-defined tensor calculus on DCq Riemannian 
manifolds, including covariant derivative of tensors, provided that the components of the tensor 
belong to BV. We restrict our attention to covariant tensors since we are mainly interested in 
defining the Hessian of a DC function on a DCq Riemannian manifold (this will be done in the 
next section). However, we will also need to consider BV vector fields as a simple instance of 
contravariant tensors in order to evaluate our tensor fields, this is done in a dedicated part. It is 
simple, on the basis of these considerations, to extend our arguments to general tensors; the details 
are left to the reader. 
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4.1. Definition of tensors on X*. 


Definition 4.1 (Covariant tensors on an open set V C M^). A covariant p-tensor S' on C with 
BV (resp. QXi) components is by definition, 

S= ^ ® ® 

ji.-" dp 

where the multi-indices {ji,--- ,jp) run into {I,--- ,N}p and S BV{V) (resp. GM{V)). 

Now, we define the pull-back of a covariant tensor with components in BV{V) or in G.M{V), 
through a transition map. 


Definition 4.2 (Pull-back of a tensor). Given F : U —> V a DCq transition map and S a 
covariant p-tensor on V with components in BV{V) (resp. t/A4(C)), we define a covariant p-tensor 
F*S on U by the formula 


(15) 


H C" )*p jl F"jp 


dFn 

dxi^ 


dF 

-^F*iSn-jp) ® ® dF^. 


Remark 4.3. F*S has BV components according to Theorem 12 .hi According to the fact that 
BVo{U) is an algebra and Remark 13.91 the components of F*S belong to BVq{U) (resp. GMoiU)) 
whenever the components of S belong to BVo{V) (resp. GMo{V)). 


Definition 4.4 (Covariant tensors on a DCq manifold X*). A covariant p-tensor S on an open 
subset fl C X* with BV (resp. GM) components is, by definition, a family 

where for each a such that UaldVt 0 , Sa is a covariant p-tensor on (j>a{Uar\Vl) with BV{(l)a{Uar\Vl)) 
(resp. t/A4((/>o,(17a nfl))) components. The family of tensors is required to satisfy the compatibility 
condition 

F*S 2 = (02 o 0 r^)* 5'2 = 5i in 01 (Ri n Rz n fl), 
as Radon measures or for points, depending on the regularity of the components, and for 

any pair of charts (Ci, 0i), (C/ 2 , 02 ) such that C/i fl C /2 fl 7 / 0 . 

We define similarly covariant p-tensors S on an open subset Gt C X* with BVo (resp. GMo) 
components. 


4.2. Covariant derivative of tensors on X*. Our next task is to define the covariant derivative 
DS of a BV tensor S on X*. We proceed as above, by defining first DS in a chart and then verifying 
the compatibility condition. 


4.2.1. Local definition of covariant derivative. To this aim, we have to introduce the Christoffel 
symbols. 


Definition 4.5 (Christoffel symbols). Let X* be a DCq manifold with Riemannian metric g. The 
Christoffel symbols are then defined as a collection of Radon measures {(ri“i)^j}agA, where 


(16) 


pfe ^ 1 ki ( d>gii dgij _ dgij \ 

2 ^ \9xj dxi dxi ) 


belongs to GM.q with singular set S, the superscript is omitted for simplicity, and gij are the 
components of the metric in the chart (Ua,4>a)- 


Remark 4.6. In the appendix, it is proved that the Christoffel symbols read in different charts 
satisfy the same compatibility relation as on a smooth Riemannian manifold. 


With this definition in hands, we can now introduce the notion of covariant derivative. 
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Definition 4.7 (Covariant derivative of a BV tensor in a chart V C R^). Given a covariant 
p-tensor S' on C = (/>(V) with BV components, we define a covariant (p + l)-tensor DS by the 
formuia 

DS = ^ ® ® dx''^, 

70 1 ■■ ■ 


where 


(17) 


{DS\o 


dS,- 


N 


dxi, 


j—1 m—1 


-pm 

ij-imijj^i ioij ' 


The components of DS beiong to QAi{V). 

In order to infer from the above definition a weii-defined notion of covariant derivative of tensor 
on any open subset of X*, it remains to prove the compatibiiity formuia 


F*{DS)=D{F*S). 


This is the goai of the next section. 


4.2.2. Covariant derivative of BV tensors on X*. We start with the foiiowing fact. 

Lemma 4.8. Let il C be an open set and f,hG BV{Ll). Assume that at least one of the two 
functions f, h belongs to yBVo(fl). Then, for any i S {1, • • • ,iV}, it holds 

djfh) ^dh ^ ^df 

dxi dxi dxi 

Proof. For mere BV{Ll) functions, the above formuia is not true in generai because of the jump 
part of the derivatives (for instance if / = /i is a characteristic function). Under the assumption of 
the iemma, at most one of the two functions has a jump part in its derivative. Then, the standard 
chain ruie formuia hoids as proved in [H Theorem 3.96] for instance. □ 

Proposition 4.9. Let F : U ^ V be a DCq transition map. Let S be a covariant p-tensor with 
components in BV{V). Then, the following equality of Radon measures holds: 

(18) F*{DS) = D{F*S). 

Proof. By iinearity of F*, and up to a permutation of the coordinates, it is sufficient to prove 
m in the case when S = Si...p dy^ 0 • • • 0 dyP. We start with F*{DS). In the foiiowing, we set 
i = (io, • • • ,ip) and j = {jo, ■ ■ ■ ,jp) for arbitrary muiti-indices in {1, • • • , By definition of 

covariant derivative. 


iDS)j= ' 


p N 


dy. 


'JO 


-EE 

k—1 m—1 


■Jk-imjk-\-i---Jp^ jojk ’ 


where the F^ stand for the Christoffei symbois in the chart reiative to V. This yieids 
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Therefore, the components of F*{DS) satisfy 


iF*iDS))^ = 


dxi„ dxi 




i 

_ dPjo ,,. 

■■■ ^^^p 




Thus, using that S = Si...p dy^ 0 • • • 0 dy^, 


(19) {F*{DS))^ = J2 


dF,„dF, 9Fp/ 


dxig dxi^ dxi. 


(-(^)) 


EE( 


dPi aPt ap^ \ aPj„ ap,^ 


dxi^ dxi^^ dxi J dxig dxi 


k=l 30,3k 

Now, we compute the components of D{F*S). Starting from 

F*S= ^■■■^S3...poFdF^ 


Si-poFF*{T%^J. 


we get. 


{d{F*S))^ 


dx 


dxi^ dxi 

.*p) 


■20 

p N 




( 20 ) (£>(^*5')) 


^)il---ik-imik+i---ip^ 

k—1 m—1 

aFi dFp d 




dxi^ dxi^ dxtg 


_vv^ oFr- 

h^^dx,, dx.,_,dx^ dx,^ '■■•" 


Now according to m in Corollary 13. Ill we have 


( 21 ) 


dxi 




Therefore, according to (I19L (l20l) . and (l2T]l . F*{DS) = D{F*S) if and only if for any mnlti-index 
i we have 


(22) -i:E( 


dFi m dFp N dFjg dF,, 


fc = l 30,jk 


dxi^ dxi,^ dxi^ ) dxig dxi^ 


\ dl<jg dtj^ o F F*(V^ i 


p N 

-EE 


dFi 

dFp\ 



dxi-, 

dxi^J 



dFi 

dFk-i 

dFk 

dFp 

dxi^ 

f)'T ■ 

ci'T 

L/U/m 

dxi^ 


S'i...„oFr 


■P " ^ ioik- 


This eqnality can be obtained from the following transformation law of the Christoffel symbols 
(whose proof in our setting is given in the appendix): 

pjTp — ^ ( 

(23) r”. = E^»P E 


ap„ ap. 


dxig dxi^ 


"■^*(rL) +E 


dFz 


■oF 


d^Ff> 


dye dxigdxi. ' 












































DC CALCULUS 


15 


Indeed, from the above formula we infer 


p N 

EE 

k—l m—1 


dFi dFk-i dFk 


dxi 


dF, 




^ s,...,OF 


p N 

= EE 

k=l 6 = 1 


dFi 

dxi^ 


dx_^ 
dFk-i m 

dxif._^ dxif. 


dFp 

dxi_ 


dFu dFy p*i-pe \ , 

dx,,dx,, ^ 


Si...poF 

d^Fe 

dxigdxi^ 


N 

(E 


-1 


dFk dF, 
Z dxm dye 


oF 


= EE 

fe = l joJk 


dFi dFk-i dFu 


dxi 


FSi...p o F 


where we use that X]m=i §i^ ° F = Ske in BVo{U) to get the first equality, and Lemma 

to get the last one. Inserting the last equality into (1^^ gives the result. □ 


dxii^_^ dxi^. 

d (dFi 


dFp dF,, dF,, 
dxi^ dxig dxi^ 


Si...poFF*{TZj 


dxig \dxi^ 


dFp 

dxi 


4.3. Vector fields. Let us recall that the notion of DCq Riemannian manifold is a particular 
instance of Lipschitz manifold. Therefore, vector fields are well-defined objects on such a manifold. 
In this part, we briefly give the definition of BV, and QAi vector fields on a DCq Riemannian 
manifold X* {BV vector fields are well-defined thanks to Theorem 12.61 and the fact that BV is an 
algebra). After that, we verify that we can evaluate QXi tensors by means of BV or Cw,o vector 
helds and get an intrinsic system of Radon measures. 


4.3.1. Cw, BV, andQAi vector fields. 

Definition 4.10 (C^, BV, and QM. vector fields). Let Vl C X* be an open set. A C™ (resp. BV, 
QAi) vector field on Vl is the datum of a system of vector fields. Namely, for any chart {U, (p) such 
that U n n 7 ^ 0, we are given a vector field 

with components in Cw{4>a{Ua H fl)) (resp. BV{(j)a{Ua H 11)), QAi{(j)a{Ua H Vl))). Moreover, we 
assume that for any pair of charts {pa, Ua), {pp. Up) such that Ua V\Up C\Vl %, the following 
compatibility condition holds 

(24) F*{Xp)-.= Y,^F*{Xl)—=Xa Af^-a.e. 

Remark 4.II. Our definition of vector field is a bit non-standard. Here, we identify vector field 
with contravariant 1-tensor. However, using that f is a biLipschitz homeomorphism differentiable 
at any point out of the range of S, it is easily seen that (El is equivalent to 



for all j G {1, • • • , N'\. In the sequel, we freely use both conditions. 

Remark 4.12 (Existence of local orthonormal frame made of Cw^o vector fields). Let (^-)i<i< 7 V 
be the constant vector fields induced by a coordinate system. These vectors obviously belong to 
Cw,o{S). By applying the Gram-Schmidt process to the (^-) with respect to g{x) at each point x 
in the domain of the chart where g is defined {i.e out of the image of 5), we get a local orthornormal 
frame {Ei) where the {Ei) are continuous out of the image of S through the chart (indeed the metric 
components satisfy this continuity property). Since the {Ei) are unit vectors with respect to g, the 
property (H on 5 yields that the {Ep are locally bounded out of a c;?^^-neghgible set. Combining 
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this with the continuity property we just mentioned shows the (Ei) are Cw,o vector fields with 
singular set the image of S. 

Lemma 4.13. Let 12 C X* be an open set and let X be a BV vector field on Ll. Then, given two 
charts Ua), Up) such that Ua H 17^ fl 12 ^ 0, if X^ is a precise representative of X read in 
the chart (</>„, Ua), then dFx{Xa{x)) is a precise representative of Xp o F. 

Proof. Since Xa coincides with its precise representative on the image of 12 \ S' through (jja and F 
is differentiable at such points, it suffices to prove that for almost every approximate jump 

point Xa of Xa, F{xa) is an approximate jump point of Xp and the one-sided approximate limits 
satisfy 

X+{F{Xa)) + X^{F{Xa)) = d,,^F{X+{Xa) + X-{Xa)). 

We shall prove this formula for approximate jump points S f2 fl Ua \ S which is sufficient for 
our purpose. 

According to Theorem 12.51 for such a Xa there exists a Borel map u : Jx^ —>■ such that 

liui'f \\Xa{x) - X^{Xa)\\ 2 dJf’^{x) =0. 

P'1'0 J Bf (Xa,v{Xa)) 

We set xp = F{xa), k{F{x)) = {*{dxF))~^{iy{x)), and v'{y) = K{y)/\\K{y)\\ 2 . We have to show 
that 

lim / \\Xp{y) - d,„F(A±(x „))||2 d.^^{y) = 0. 

PfO JBfix^yixp)) 

To this aim, we point out that F~^{Bp{xp)) C BBip(F-^)p{Xa) and F~^{y) G B^{xa,v{xa)) 
whenever y G B^{xp, v'(xp)) and p is sufficiently small. Indeed, since xp is a regular point by 
assumption, F~^ is differentiable at xp; hence 


{p- {y)-F- {xp),v{xa)) = {dxpP- {jj-xp),u[xa))+o{\\y-xp\\ 2 ) 

= {y- xp,* {dxpF-^)u{Xa)) + o{\\y - xp\\ 2 ) 
= \\Kixp)\\ 2 {y - xp,v'{xp)) + o{\\y - xph). 
Therefore, applying the change of variable formula, we get 




Bf {xpF'i^p)) 


\\Xp{y)-dx„^F{X^{xa))\\2d^^{y) < 


L 






< (LipF)^ [ 

J E 


\\dxF{Xa{,x)) - dx^F{X^{Xa))\\2d^^ {X) 


by . "We conclude by using that dxF is continuous at Xa and Xa i® an approximate jump point 
of Xa- □ 

4.3.2. Covariant derivative ofBV vector fields. As for covariant tensors, we start with the definition 
when read in a chart. 

Definition 4.14 (Covariant derivative of a BV vector field in a chart V C K^). Given a vector 
field Y on V = 4){V) with BV components, we define a (1, l)-tensor DY by the formula 

DY = Y.iDY)]dy^O^^, 

j,S 
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where 

(26) 


dY, 


N 


11=1 


The components of DY belong to QM.{V). 

Definition 4.15 (Pull-back of (1, l)-tensor). Given a DCq transition map F : U —> V and a 
(1, l)-tensor 5" a on 1^ 


with components Sj in BV{V) (resp. QA4{V)), we define a (1, l)-tensor i^*5 on U by the formula 
(27) 


F*S= V ^^^oFF*{S;)dF^^. 

oxi ays oxk 


We refer to m for a general definition for (p, fc)-tensors. As for covariant tensors, in order to 
infer from this local definition a well-defined notion of covariant derivative of BV vector field on 
any open subset of X*, it remains to prove a compatibility formula. 

Proposition 4.16. Let F : U ^ V be a DCq transition map. Let Y be a vector field with 
components in BV{V). Then, the following compatibility formula holds 

(28) D{F*Y) = F*{DY). 


The proof is very similar to that of ProDOsition l4.9l thus the details are left to the reader. 


4.4. Norm of tensors. In order to define the norm of tensor with Radon measures components, 
we cannot apply the standard definition for smooth tensors on smooth Riemannian manifolds (see 
for instance m Section 6.3]). Indeed, the classical definition involves products of components of 
the tensor and thus does not apply to our setting. However, in the case when the tensor components 
are functions, our definition coincides with the standard one. 


4.4.1. Evaluation of tensors. 

Definition 4.17 (Evaluation of tensors). Given H C X* open, a covariant p-tensor S with QAi 
components on H, and Xi, ■ ■ ■ ,Xp vector fields with Cw,oi^) components, we may evaluate the 
tensor on the vector fields by the formula in any local coordinates 

sp{xpp, • • •, Xp^p) = ^p,i ■ ■ ■ 
i 

Note that it is also possible to evaluate tensors by BV vector fields. However, it is necessary to 
use the precise representatives of the vector fields for the formula to make sense, since Spj need 
not be absolutely continuous w.r.t. . Whatever the regularity of the vector fields, what really 
matters is that they are defined everywhere out of a J^^“^-negligible set and bounded in the sense 
of Definition 12.91 

The following lemma shows that this definition satisfies a compatibility condition. 

Lemma 4.18. Suppose we are given a covariant p-tensor S with QA4 components defined on an 
open set D C X* and BVor C^.o vector fields Xi, ■ ■ ■ ,Xp on D. Then, for any pair of charts 
Ua), ((/>/ 3 , Up) such that Uct n n D ^ 0 , the following compatibility conditions holds 

F*iSpiXpp,--- ,Xp^p)) = s4x^p,--- ,x^,p) mf^^iu^nupnn). 

Therefore the collection of measures {Sa{Xap, ■ ■ ■ , Xa^p))aeA in GA4 is a system of Radon mea¬ 
sures on LI according to Definition \S.l‘A 
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Remark 4.19. Note that given / a BV function and X a BV vector field defined on 17, the above 
result and Lemma [11.141 yield that df{X) is a well-defined Radon measure on 17. 


Proof. By definition, 


SpiXp,,, • • • , X0,p) = Y, 


'Jp ’ 


In other terms, Sp{Xp^i, ■ ■ ■ ,Xp^p) is a sum of measures where each summand is a product of a 
locally bounded Borel function and a Radon measure in QM.{(j)a{UaP\UpPiVl)). Therefore, applying 
m, we infer 

F*{SpiXp,,, ■ ■ ■ ° ^) ■ ■ ■ ° F)F*{SpY 

i 

Now, by definition of vector field (see (I^Kll l and according to Lemma [4.131 we have the following 
equality up to a Jf'^“^-negligible set 

E 


Besides, by definition of F* on tensors, 

dF 
dx 


_ firr. a,k- 




By combining all of this, we hnally get 

• • • , Xp,,)) = Y^P*^Sp)\_Xl, ■■■XI, 


a.,p 


i 

— Bol {^a,l ; * * * ) ^a,p); 


and the proof is complete. 


□ 


Remark 4.20. It is also possible to evaluate a covariant 1-tensor with BV components on a vector 
field with QAi components. The above proof applies with minor changes. 

4.4.2. Local definition of the norm of tensor. We first define the norm of a tensor in a chart. 
Recall that the existence of local orthonormal frames made of Cw,o{S) vector fields is guaranteed 
by Remark 14.121 

Definition 4.21 (Local definition of |S'|g). Let 5* be a covariant p-tensor defined on an open 
set 17 C X* with SA4(fl) components and {Ua,4>a) a chart of X*. Let (£'i)i<i<Ar be a local 
orthonormal frame in (f)a{Ll fl Ua) with Ei G Cw,o{fia{^ H t/a), fiaiS H !/«))■ We define the norm 
\Sa\g of Sa as the total variation of the -valued Radon measure defined as 

(29) ('S'(i?ii, • • • ,i?ip))(ii,...,jp)g{q...,Ar}P. 

Here and in the sequel, we always assume that local orthonormal frames are C^.o regular out 
of the image of S through the chart. Next, we show that Definition 14.211 is well posed. 

Lemma 4.22. Let S be a covariant p-tensor defined on 17 C X* with QM{Ll) components and let 
(Uonfia) be a chart of X*. Then, the norm liSctlg does not depend on the choice of the local frame. 

Proof. Let (i?i)i<i<Ar, (i?i)i<i<A be local orthonormal frames on 17 and let P : 17 —>■ be the 
change-of-basis matrix with respect to {Efi and {Ei), i.e. Ef = J^eFuEj. We denote by ct = |A| 
the total variation of the vector-valued measure in (EH) and by a and A the analogous quantities 
for Ei. Our goal is to show that a and a coincide as measures in 17. To this aim, we will find a 
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representation of A in terms of A. Let T : Ai —>■ be a Borel unit vector field providing the polar 

representation of A, i.e. 

,£'jp) =T^(t for alli= (ii,--- ,*p) G {I,-- - ,N}p 

and, writing S = Sj (8) • • • 0 , let us compute: 

i 

= E 4' 

N 

= X! ■ ■ ■ Pipip^iEe^ , • • • ,Ei^) 

h,--- ,£p=i 

N 

= X! Pilii-■ ■ PipipTii,--,ep<^- 

,£p=i 

Setting {P*)u = P^ii ■ ■ ■ Pt^ip, this proves that 

\ = Y,iP*WL(T- 

I 

Since P is an orthogonal matrix and ||T(a ;)||2 = 1, it is easily seen that ||P*r(x )||2 = 1, therefore 

|A|=ct. □ 

We can now check that is a well-defined Radon measure on Al. 

Lemma 4.23. Let S be a covariant p-tensor defined on Cl C. X* with QM{Cl) components. Then 
(l^als) ctgA is a system of QM measures. Thus, it induces a measure G QM{Cl) defined by 
called the norm of S. 

Proof. By definition of tensor, for any pair of charts {(pa, Ua), {Pp, Up) such that UaL\UpC\Ct %, 
the following compatibility conditions holds 

F*{Sp{Ep^l,--- ,Ep^j,)) = Sa{Ea,l,--- ,Ea,p) 

where {Eap, ■ ■ ■ , Ea,N) and {Ep^i, • • • , Ep^f^) are local orthormal frames in pa{UaL\Cl) and pp{UpC\ 
n) respectively. The result then follows from Lemmas 13.131 and 13.141 □ 

Remark 4.24. Let / G BV{Ct), being Cl C X* an open set. Then, the differential of / is a covariant 
1-tensor with QXi{Cl) components. As a consequence, the above result applies and gives us an 
intrinsic notion of “total variation” \df\g of the measure-valued 1-form df. 

The arguments above can easily be adapted to show that there is a well-defined notion of norm 
for vector field with QA4 components. We give the formal definition below and leave the details to 
the reader. 

Definition 4.25 (Norm of a vector field with QA4 components). Let A be a vector field with 
QAi{Cl) components, being Cl C X* an open set. Then the norm \X\g G QfA{Cl) of X is defined 
as the nonnegative measure induced by the following system of measures. These local measures 
are defined in each chart whose domain intersects Cl as the total variation of the M^-valued Radon 
measure {6^,--- ,0^) where (0*)i<i<Ar are the local coordinates of X with respect to a local 
orthonormal frame. This definition does not depend on the choice of the orthonormal frame. 

As in the smooth setting, the norm of tensor can be used in the following estimate. 
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Lemma 4.26. Suppose we are given a covariant p-tensor S with QM components defined on an 
open set n C X* and BV or 0^,0 vector fields Xi, • • • , Xp on it. Then, the following inequality 
between measures holds 

\S{X,,--- ,Xp)\<\X,\g---\Xp\g\S\g, 

where \Xi\g = ^/g{X^,Xi). 


Proof. Since we proved that ,Xp)\ and are both systems of Radon measures, we 

can work in local coordinates {U, ff) omitting the -a for simplicity. Let {Ei) be a local orthonormal 
frame. Then, according to Definition 14.211 the Radon-Nikodym theorem implies the existence of a 
]R^*’-valued map M{x) with unit Euclidean norm such that 


Decomposing Xi, - ■ ■ , Xp in the basis (Ei) as Xg = yields for any Borel set A 

,Xp)(Al)| = | 0\----e;^S{E,,,--- ,E,^){A)\ 


■ ■ ■ I'ip 


< 


< 



■■■6l^{x)M,,...,.^{x)d\S\g{x) 

■■■\ei^{x)\\M,,...,,^\{x)d\S\g{x) 


^ Ig ■ ■ ■ 

J A 

<|Xi|,..-|Xp|,|,S|g(^), 


where to get the penultimate inequality we use the Cauchy-Schwarz inequality, the fact that M{x) 
has unit Euclidean norm, and that {Ei) is a local orthormal frame for g, hence = l-^slg- 

Consequently, since A is arbitrary, we get the result by the definition of total variation. □ 


5. Hessian and Laplacian of a EC function 


In this section, we use the tensor calculus we introduced in order to define the Hessian, its norm, 
and the Laplacian of a EC function. Let / : H C X* —R (with D an open set) be a EC function. 
Then df is a covariant 1-tensor with BV components whose representation in a chart ([/, 4>) is 


d{fo(j) 1) 




5.1. Hessian of a EC function. We can define the Hessian of / as follows. 


Definition 5.1 (Hessian of a EC function). Let X* be a ECq Riemannian manifold and let / 
be a EC function defined on an open subset of X*. We define the Hessian of / as the covariant 
2-tensor 

Hess/ = Edf. 

When read in a chart {U,(j)), the components of Hess/ are given by 


(Hess0/)ij 


^) _ difocf ^) fc 
dxidxj ^ dxk 

k—1 


The norm of Hess/ is then defined as in Lemma 14.231 


Remark 5.2. Note that the components (Hess^/)^ of Hess/ in local coordinates are symmetric 
with respect to i, j. Therefore, the Hessian of a EC function can also be considered as a symmetric 
matrix-valued Radon measure. 
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Proposition 5.3. Let Q C X* be an open set, f € DC(fl), and X,Yg Then, the 

following equality in QM(Ll) holds 

(30) Hessf{X,Y) = D{df{Y)){X) - df{DxY). 

Remark 5.4. df is a covariant 1-tensor with BV components, therefore it makes sense to evaluate 
it on a QAi vector field, see Remark |4. 201 D{df(Y)) stands for the derivative of the BV function 
df{Y) (which is usually not written d{df{Y)) to avoid confusion with the exterior derivative of 
differential forms). 


Proof. It is sufficient to check that both measures coincide when read in a chart. For simplicity, 
we keep the same notations for the function and the vector fields read in a chart. Let us recall 
that read in a chart, 


(31) Hess/(X,y)= ^ 

while the right-hand side of ([5n|) is the sum of 


dxidxj 


N Of 

^dxk 




(32) 


D{dfiY))iX)= 


dxi\dxn 


and, using DxY = Ei<*.j<iv ^*(S + A 

(33) 




From (ED, ED and ESD, it is then clear that the absolutely continuous part and the Cantor 
part of both measures coincide since 

D'^{uv) = uD^v + vD''u 

when u, v are BV functions and stands for either the absolutely continuous part or the Cantor 
part of the derivative, we refer to [4] for a proof. Therefore, it remains to study the jump part of 
the measures. Let us recall that, according to Theorem 12.5l dl. the jump part of the derivative of 
f € BV has the following structure 

D^'^f = {f+-f-)ujdJY’^-^LJf 

where Uf : Jf ^ is a Borel function, the so-called approximate unit normal, and Jf is the 

set of approximate jump points of /. 

To this aim, we fix i, j G {!,••• , N} and compute the jump part of {^^Yj) appearing in 

(ED). According to the chain rule formula, we have to consider three cases. Below, we write -1“ 
for the jump part of the partial derivative and (•, •) for the standard Euclidean inner product. We 
also use many times that we are considering the precise representatives of Yj and which are 
“undefined” only on the J^^“^-negligible set of approximate discontinuity points which are not 
approximate jump points. 

First, on \ Jgf/Qxj one has 


= §^iY;-Y-){d/dx.,l.Y,)dJff’^-^ViJy^\Jof/0.,). 


dxi ^dxi 


Second, on Jg/jdxj \ Jvj one has 


'(.o.r.. y) dxj ^ Idxj) dJif L {Jdf/dxj \dYj)- 


df 


dxi ^dxi 


Last, assuming with no loss of generality that uy. = lygf/dxj ^-a.e. on Jgf/dxj C Jy (here 

we use the fact that the approximate unit normals coincide up to the sign on the 
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intersection of approximate jump points, see [U Example 3.97], and so we can assume that they 
coincide, up to a permutation of the right and left approximate limits), on Jojjdx^ H one has 


Qju 

dxi 


(Ky) = (Y+^ 


dxi 




) {djdxi.vQf/Q^.) ^\_{jQf/Q^. (IJy). 


According to (15^ . ((M)) . we have to compute the density 9 with respect to ^ of 


^(Ky-)-^^Y 

dxidxj dxj dxi ^ 

On Jy. \ Jdf/dxj we have 

" = - y-) - - rp) 0/3-.- = »■ 

On Jdf/dxj \ Jvj we have 


Finally, on Jgf/Qxj O Jy^, writing 0 = 9{djdxi^VQf = 0{djdxi,VYj), on this intersection we 
have 


s-_USL\,.iL-r-Yiiy;-ypilL\IL 

df 


V ^ dxj ^ dxj ^ 


(y.+ + Y-) - (Y+ + Yr 

2\dxj dxj ^ ^ ^ 


V'^5a;, dx, ^ V' 


Therefore in each case the measure 


^(^Y)-^ — Y- 
dxidxj dxj dxi ^ 


r)^ f 

coincides with the jump part of Yj ^ ./ . ■ Since i, j are arbitrary, the proof is complete. 


□ 


5.2. Laplacian of a DC function. In order to define A®/, we proceed as for tensors by defining 
the Laplacian locally in a chart first and then verifying the required compatibility formula. 


Definition 5.5 (Local definition of A®/). Let / : fl C X* — >■ M be a DC function with O an open 
set. We define the Laplacian A^/ of / read in the chart (p as the trace of tless^f with respect 
to g. More precisely, if {U,(j)) is a chart of A* such that t/ fl fl ^ 0, in the system of coordinates 
induced by (p one defines 


A^/:= ^ 


I dxidxj 


^d{fo^ \ 

k V 


in p{u n n). 


Now, we prove that the local definition A^/ of A®/ provided above induces a system of Radon 
measures. 


Lemma 5.6. Let F : U —> V be a DCq transition map relative to the charts {U,p) and (V,^). 
Let f : n — > K &e a DC function and suppose ff fl R fl R ^0. Then, 


F*{A^f) = A^f. 


( 34 ) 
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Proof. Note that in the following argument, the Hessian is considered as symmetric matrix-valued 
Radon measure. By dehnition, 

F*(A^/)=F*(tr(G-iHess^/)). 

We set F*((Hess^/)..) the symmetric matrix-valued Radon measure whose entries are R*((Hess^/)ij) 
in the sense of Definition 13.81 With this notation, we can rephrase the above equality in the fol¬ 
lowing way 

(35) F*{A^f) = tr (g-1 o FR*((Hess^/)..)) . 

From the tensor equality Hess^/ = F*(Hess^/), we infer the following equality of matrices 

Hess^/ =‘ (dF)F*((Hessj/)..)(dF), 

which is equivalent to 

R*((Hess^/)..) = \dF)-^mss^f{dF)-\ 

By combining this together with 

G-^oF = {dF) (G-^YidF) 

and then inserting these equalities in (I35L we obtain 

F*{A^f) = tr{{G-^)Eess^f) 

and the proof is complete. □ 


Now, using Lemma l5.fil and Lemma |3.141 we can define the Laplacian of a DG function on D. 

Definition 5.7 (Laplacian). Let / : D C X* — ?> R be a DC function. We set A®/ the Radon 
measure defined by m- 

Our goal is now to prove that A® / coincides with the weak Laplacian defined through integration 
by parts. To this aim we generalize to our setting the following classical formula. 


Proposition 5.8. Let X* be a DCq Riemannian manifold and f be a DC function defined on an 
open subset of X*. Then the local expression A^f o/A®/ in a chart (U^cj)) is given by 

1 _ a / _ 

(36) A^f = 


Vdet G 




Proof. Let us set f = f o (p ^ for notational simplicity. We expand 








I dga 

\ dxi dxj dxi 


2 dx, 


-9^9 


7 

1 dgij 

\dxi ) 

2 dxi _ 


= E 

since g'^^ = g^'". By differentiating with respect to Xi the equality 'Ylii<j<N 9^^^93i ~ '®^ 0 j we 

infer 






V dxi 


E 


l<i,j,k,l<N "■ \ / l<i,j,k,l<N 

so that the definition of A^f gives 

ij J 


dxJ ^ 


dxj dxi 


(37) 


A^/ = ^ F- 






2 dxk ^ dxi 


E 
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Now, we expand the right-hand term in (I5S1) using Lemma IT51 and then apply ([S]) to get 


VdiFC 


E 


l<i,j<N 


A 

dxi 




l<iJ<N 


d^f 

dxidxj 


E 

l<i,j<N 


df dg'-^ 

dxi dxi 


E ^ 


df 1 aVdiFG 

dxj VdetG dxi 
df 


^ dxidxj ^ dxj dxi 


df 

r\ y y r\ 

OXj OXj 

l<i,j^k,s<N ^ 


E 


and the proof can be completed comparing with (13711 . □ 

Proposition 5.9. Let f : LI C X* —> R be a DC function and ip G Lipc{Lt). Then, 

[ g{yf,yip)dvg = - [ ipA^f. 

J fl J Q 

Proof. By using a Lipschitz and locally finite partition of unity = 1 subordinate to the 

atlas of X*, we have 


[ i/'A9/ = V f ipo y/ det G A^/ 


S / , ° ^ ^ AP) 


= -T,f 

= -E / dc.g{X^P.Xf)dvg + 0 

aGA-^f^ 

where, to get the last equality, we use the fact that the partition is locallly finite and thus for 
Jf'^-a.e. xGLI, V(X„eA Oa){.x) =0. □ 

6. Integration by parts formula for the Hessian 

In this section, we prove that the Hessian of a DC function satisfies the same integration by 
part formula as in the r 2 calculus, see [6] for a precise definition. 

To this aim, we generalize to our setting the following classical formula. 

Proposition 6.1. Let ip be a DC function. The following equality of Radon measures holds 

D^^ViP=^V{\V^\l) 

Proof. It suffices to prove the result when the objects are read in a chart {U,(p). However for 
simplicity, we drop o and simply write / and ip in the computations below. Let us recall that 
in a chart, the covariant derivative of a vector field Y with BV components is defined as 

DY = Y,iDY),,kdx^ - ^ 

i.k 


dxk 


where 

(38) 


N 
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Decomposing Vip in the standard basis, we get 

VV’ = 3 

This leads to 

g 


l<k,s<N 


ks 


ks dj’ d 
dxs dxk 


l<k,s<N 


D 3 V'd’ 

OXf, 


E 


l<k,s,l<N 


dip 

dx/. 


d 

dxk 

( S “'"S) 


\l<t<N ^ j 


E 5 

l<t<N 


dxf ^ 


Now, if f, g G BV{Vl) are such that at least one of the two functions f,g belongs to ,BVo(D), 
then for any i G {1, • • • , N}, it holds 


This yields 
(39) 
where 


dxi dxi ^ dxi' 


= A + B, 

dp: d'^ip d 


E 

l<k,s,l,t<N 


gkSglt 


B= E 


dxs dxkdxt dxi 
ks dp) dp: (dg^* d 


l<k,s,l,t<N 
d 


dxs dxt V dxk dxi 


+ g^*D a ) . 

^ ^9X1J 


Now, we compute for a fixed i £ {1, • • • , A^}. We notice that 


W „fc5 9p: dp: ^ 
9\B) g^.) — / ^ 9 


l<k,s,l,t<N 


dxs dxt 


where 

Using that 
(40) 


Ci,t,i = ^9H+9'* E r 


dx 


kl 9ui ■ 


l<u<N 


9u9‘* = <5*1 


l<i<A 


hence a constant function (here S stands for the Kronecker symbol) whose derivative is null, we 
can rewrite 

dgii 


kl 9u 


E E E n 

1<Z<7V 1<1<N ^ l<u<N 

Now, since T'^i — ^ + ^ - ^))i get by expanding X) 1<U<N ^kl 9‘U 

^ }_ It ( ^9ki dgii dgki \ 

2_^ i,t,i 2_^ 2 ^ V dxi dxk dxi) 


KKN 


KKN 
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We filially obtain 


E 

l<k,s,l,t<N 


E 

l<k,s,l,t<N 


1 ks It dj’ d-ip /dgki 

2 ^ dxs dxt V dxi 

ks It dj’ dgu 

2 ^ ^ dxs dxt dxi 


dgii 

dxk 


d9ki \ 
dxi ) 


since the other terms cancel because of the symmetries. Using (001) again yields 


(41) 




E 

l<t,s<Af 


1 dg*^ 

2 dxs dxt dxi 


The next step is to compute 


1 aivy-ls 

2 dxi 


Starting from the equality 


iViAl^ 


E 

l<t,s<N 


dip dtp 
dxs dxt ’ 


we infer 
(42) 


1 

2 axi 



By comparing with ((HU and ((IT]) , we get 


(43) giD^^VtP, ^) - i ^ = Ei<....<a 5 


^ks dj) 

dxs dxkdxi 


_ 1 ^fcs d ( d'lp dtp \ 

2 l^l<k,s<N y dxi \dxs dxkj 


where we use (iTOI) to simplify the first term on the right-hand side. Now, if tp were a DCq function 
then it would be no jump part in its second derivatives and we would get 


d / dtp dtp \ 

dxi\dxsdxk^ dxidxs dxk ' dxidxk dxg 


d^tp dtp ^ d^tp dtp 


Since the second distributional derivative of a function is a symmetric matrix-valued Radon mea¬ 
sure, the result is proved in this special case. But tp is not DCq in general, thus we have to 
compute the jump part of the term in the right-hand side of (1431) and prove that it vanishes. By 
symmetry of the second derivative, for any t, s S {1, • • • , N}, the jump part of and 

coincide. Let us also recall properties of the jump part ; in the following W^tp stands for the 
Euclidean gradient of tp. First, for point in J sv- C\ J , v = v di, = uye,* and 

idxp, idxs i^Xg 

up to a negligible set J . In what follows, we make the convention that 

= 4^ = 4^ out of J , SO that the precise representative satisfies 

^ fc ^ k ^ k y a; T, 


dtp 1 / dtp ^ dtp \ 
dxk 2\dxk dxk ' 

everywhere out of a negligible set. This convention allows us to write for any k £ 

{I,--- ,iV}, 


/ dtp dtp 

V dxk dxk 




^ dip 
\ dxk 


S') K"*, 


'U, 














DC CALCULUS 


27 


Using these properties, we can now compare the jump parts of the derivatives. Starting with 


E ^ 

l<k,s,<N 


ks 


d'lp 

dxs dxkdxi 


= E ^ 

l<k,s,<N 

= E ^ 

l<k,s,<N 

- ( i: 

l<k,s,<N 


. dijj 

dxs dxidxk 


r dip ^ dip \ 

/ dip + 

dj’ \ 

\dxs dxs J 

\dxk 

dxk ) 


0 '=* / dijj 


+ 


dip 

dx. 


+ 

dip \\ 

\dxk 

dXk 7/ 




- ( 

g'^‘1 ( dip 

X 2 (Z 



l<k,s,<N 

We now consider 


l<k,s<N 

ks d ( dip dip \ 
^ dxi dxs dxk j ' 

On J dtp 0 J dip 5 

M 

dip d'lp 
dxs dxk 

\ _ f dip dip 

) ~ dxk 

On J dip \ J dip , 

dxk dxs 

dxi \ 

dip dip 
dxs dxk 

+ 

II 




_ \(d±+ , _ 
2 dxs ' ( 

On J dip \ J dip , 

dxg 

dxi \ 

dip dip 
dxs dxk ^ 

\ _ dip f dip 

/ “ dxk y dxs 

- if _L . 

2 1 dxk '' 1 


+ dip ~ dip' 
: dXs dXk 


\{yp^a)^^ ^Lj 






d'lp 

dxk 


d'lp 

dxs 


)0 


’ - / ’ ' 




■)(^ 


dip 

dxk 


dip dip 

dxs dxs 


dxk 

^ Lj, 

Consequently, we get thanks to the symmetry with respect to k and s, 

1 V gfc^ / + dip~ dip~\\ 

^ ^ dxj V dx. dxk / V ^ 2 V dx. dxt dx. dxk )) 


V^i{j 




^V^'lp, 


a )Jv^ 


l<k,s<N 

The equality is proved. 


l<k,s,<N 


□ 


With the above result at our disposal, we can now establish the integration by parts formula 
involving the Hessian. 

Proposition 6.2 (Integration by parts formula). Let v (respectively u) be a DC (resp. DCq) 
function defined on an open subset H of X*. Then, for any compactly supported Lipschitz function 
Ip defined on LI, the following equality holds 

I Ip Hessv{Vu,Vu) = — I ipg{Vv,Vu)A^u—l- j ipg{Vv,X\Vu\‘^)— j g{Vv,Xu)g{Xu,Vip) dvg. 
Jq Jq, ^ Jq Jq 

Proof. The claim follows from the integration by parts formula 

( 44 ) / Ip g{Vv,Vu)A^u = — / g{Vu,V{g{Vv,Vu)ip)) 

Jn Jq 

proved below. Let us first explain how to complete the proof from this equality. First, notice that 
g(fi7v,Xu) and ip are BV functions on Ll and that ip has no jump part in its derivative. Therefore 
according to Lemma H751 the Leibniz rule holds for these functions: 


X{g{Xv, Vu)ip) = ipV{g{Vv, Xu)) + g(Xv, Xu)XiPjr‘^. 


This yields 


( 45 ) 


/ ipg{yv,Vu)A^u = - [ ipg{yu,X{g{Xv,Xu))- f g(yv,Vu)g{Vu,Xip)djr^. 

J n J n «r2 
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Now we can rewrite the first term on the right-hand side as 


[ tljg{yu,V{g{\7v,Vu)) 
Jq 


(46) 


[ tp D{dv{yu)){Vu) 

Jq 

/ i/iHessi;(Vu, Vu) - 

J Q, J Q 

f '0Hessi;(Vu, Vu)-I-- [ tp g{yv,V\Vu\1) 
Jq 2 Jq 


p) g{Vv,D^u'^u). 


where we used Proposition 15.31 to get the second equality, and Proposition 16.II to get the last one. 
Inserting (1461) into ((151) gives the result. 

It remains to prove (HU) . Reasoning as in the proof of Proposition 5.9, we can further assume 
that ip is supported in the domain of a chart, thus it suffices to prove the result in local coordinates 
defined on an open subset of K.^ that we also called If for simplicity. Let (p£)£>o be a family 
of standard radial mollifiers. Let us set h = (/(V?;,Vu) and note that h G C'u,(n). Therefore, 
h * pg —>■ h pointwise out of a cr-finite set with respect to The Lebesgue dominated 

convergence theorem then yields 


(47) 


/ ^g(Vw, Vm)A®m = lim / tph*psA^u 
Jq efo Jq 

= lim ( / Ip g{\7h * pi;,\7u)dvg + / h * p^ g{\7'tp,\7u)dv 

AO \Jq Jq 


Vn. 


Since \g{'\/^p,'Vu)\ is compactly supported and bounded, it is clear that 

lim / h* psg{yip,Vu)dvg = / giV v u) giV ip u)d' 

AO Jq Jq 

We prove the convergence of the remaining integral thanks to Proposition Id.lOl Let us check that 

the hypotheses of Proposition 13.101 are satisfied. Using the coordinate system (xi)i<i<N^ we get 

g{Vh* ps,Vu) = 


2=1 


dxj 


By assumption, g’’^-^ip G C'u,,o(f2) while ^ G Moreover the convergence 

in the duality with C'c(U) and (IIOI) hold (see for instance [H Theorem 2.2]). Therefore 

lim / i/) g(Vh * pe, Vm) = / p!g(Vh,'S/u) 

AO Jq Jq 

and the proof of (Hill is complete. 


dh 

dxj 


□ 


7. Alexandrov spaces 

In this part we prove that our results apply to an open dense subset A* of a finite dimensional 
Alexandrov space (A, d) with curvature bounded from below. For an introduction to these spaces 
we refer to the book [7], see also [2] where all the notions below are discussed with further details. 

7.1. The quasiregular set of an Alexandrov space is a DC Riemannian manifold. 

Definition 7.1 (Quasiregular set A*). Let (A,d) be an A-dimensional Alexandrov space with 
curvature bounded below by k. Given S > 0, a point x is d-regular if there exist A pairs of points 
(pi,gi),--- ,ipN,qN) such that 

{ Zpixpj > ^ — S for alH ^ j, 

ZpiXQi > TT — S for all i 

where /.xyz stands for the angle at p of a comparison triangle in the space form of curvature k. 
The collection of pairs {pi-,qi) ■ ■ ■ (pN^qN) is called a d-strainer (at x). By lower semicontinuity 
of angles with respect to x, the set of d-regular points is open. We shall denote by A* the set 
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of quasiregular points, namely the set of (iAr-regular points for a fixed 6n ^ 1/-^, omitting the 
dependence oi Sn for simplicity of notation. 

Remark 7.2. It is proved in [7] (see also the original paper El) that 5n = iQQj^ is a suitable 
choice. Namely, for such a choice, the open set X* is a dense subset of X and a Lipschitz manifold. 
Note that it can be proved that these properties remain true for the set of i5-regular points where 
0 < S < Sjv (see [71 Corollary 10.8.24]). 

Definition 7.3 {Reg{X) and Sing{X)). The set Reg{X) of regular points is the set of points 
whose tangent cone is isometric to dimensional Euclidean space. Equivalently, it is the set of 
points which are <5-regular for any i5 > 0. The complement in X of Reg{X) is called singular set 
and denoted by Sing{X). 

Later, it was proved that actually X* is a DCq Riemannian manifold according to our termi¬ 
nology. 

Theorem 7.4. Given {X,d) a N-dimensional Alexandrov space with curvature bounded from be¬ 
low, its quasiregular set X* is a DCq Riemannian manifold with singular set S = X* fl Sing{X). 
More precisely, there exist a Riemannian metric g and locally biLipschitz charts (j) : —>■ M" 

defined by the formula 

4’ ~ i^Pl ) ■ ■ ■ ) '^Pn ) 

where 

dpiix) = -f d^,{x)dM’^{zi), 

J B{pi,ei) 

such that X* is a DCq Riemannian manifold. Moreover, the components gij of the Riemannian 
metric when read in a chart, belong to BVo and satisfy 

(48) ^Iblli < - ^11^112 forallp€R^,forx€^{U'\S), 

id 

with c = c(U') > 0 for all U' (s U. Last, the components gij are differentiable Lebesgue almost 
everywhere. 

Proof. Property (a) of the definition of DCq Riemannian manifold was proved by Otsu and Shioya 
[TB] (with the exception of the DCq character of the transition maps which is due to Perelman 
|17| 1 as well as the existence of a Riemannian metric g whose components gij are in out of 
S = Sing{X). The BV character as well as the differentiability property of gij was proved by 
Perelman in HZ]. The estimate (I^Rll is a consequence of results proved in these papers but is not 
properly stated as such. Let us give a proof of it for the sake of completeness. In the following 
(•, ■)x stands for the inner product defined on the tangent cone (note that, strictly speaking, it 
is an inner product only when x £ Reg{X) and in this case it coincides with g at x) by the formula 

(u, v)x = |u| |'c| cos Lu, V. 

Now, any fi := dp. is a semiconcave function, thus it admits directional derivative (see [19j l 
along any tangent vector u £ K^. Moreover, the directional derivative satisfies for all x £ Lj and 
u € Kx 

fl{x,u) < {fi[x),u)x 

where fi{x) is the “gradient” of fi at x (see m for more details). Using the definition of fi as an 
average of distance functions, one can prove that 

for all M £ RTa;, f'{x,u) = {^i{x),u)x 

whenever x £ Reg{X). 

Indeed, given a unitary geodesic 7 starting at x, the first variation formula (and Lebesgue’s 
dominated convergence theorem) yields 

Ml{s)) = Mx) Psi cosZ7'(0),tS dJf^ip) -f Ox{s) 

J B(pi,ei) 
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where fS stands for the unit direction of the (unique for Jf'^-a.e. p S B{pi,ei)) geodesic from x 
to p. See m for more details. Consequently, 

/ il'iO)Ax)xdJf^{p) = i cosZ 7 '( 0 ),tS dJf^ip) = /■(x, 7 '( 0 )) < ( 6 (a:), 7 '( 0 ))x 

J B{pi,€i) J B{pi,ei) 

and, by density of the geodesic directions 7'(0) in the space of directions Yjx{X) (the unit sphere 
of Kx), the above formula holds for any u G K^- Since the extreme terms are both linear with 
respect to u when x G Reg{X), we do have equality in the inequality above whenever x G Reg{X). 

Let now U' (s U. Using this equality for all i G {1, • • • ,X} and the fact that 0 is a locally 
biLipschitz map, namely the existence of c = c(l7') > 0 satisfying 

^d(x,g) < ||^(a;) - d{y)\\2 < cd{x,y) Vx, y G U\ 

we get, for all x G Reg{X) DU', u G 

(49) - < ||((Ci(a;),u)a;)l<i<A||2 < C, 

c 


since 4>'{x,u) = {{^i{x),u)x)i<i<N- We also infer from the above formula that {^i{x))i<i<N forms 
a basis of whenever a; is a regular point. In particular for all p G we can find Up G such 
that 

(j)'{x,Up) =p. 


Then, (HHl) yields 

-\J{Up,Up)x < ||p||2 < C 
by homogeneity and the proof is complete. 


(Up, Up)x 


□ 


7.2. Almost everywhere second order Taylor expansion for DC functions. In this part, 
following Perelman’s ideas CZ], we fully extend to the setting of Alexandrov spaces the classical 
result of Alexandrov on the existence of second order Taylor expansion for a convex function at 
Lebesgue almost every point. To be more specific, our goal is to prove the following result. 

Proposition 7.5. Let f be a DC function defined on Lt C X*. Then, for -a.e. point x, 


f{y) = f{x) + \xy\df{x){-tl) + i|a:?/pi7ess“‘=/(t^,tD 

where \xy\ = d(x,y) and denotes the direction at x of an arbitrary geodesic from x to y. 

Remark 7.6. We don’t know whether this result holds on any DCq Riemannian manifold. The 
proof we give below follows the arguments given by Perelman in dZl Proposition p7] with more 
details added. Indeed, this result is of local nature, so there is no need to have a well-defined global 
object Hess to consider this question. However, as stated, it shows that our notion of Hess gives 
geometric informations on the function. 

To prove the result above, we need first to introduce normal coordinates. 

7.2.1. Normal coordinates. In this part, following Perelman’s ideas, we establish some properties 
of normal coordinates around “good” points of Alexandrov space. We think these results could be 
useful in other contexts, so we state them independently of Proposition 17.51 

Remark 7.7. Note that by definition of a DCq Riemannian manifold, if the metric components are 
differentiable at point F(x), x G Reg{X), when read in a chart F then this is also true in any other 
chart whose domain contains x. The same property holds for the second order Taylor expansion 
of any DC function defined on an open subset of A*. Consequently, it makes sense to say that the 
metric g (or a mere function) is differentiable at a point x G Reg{X). According to Theorem 17.41 
g is differentiable at .^'^-almost every point. 
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Proposition 7.8. Consider a point x G Reg{X) where the metric is differentiable. Then, there 
exists a normal coordinate system N defined on a neighborhood of x and compatible with the DCq 
structure (i.e. it makes a DCq mapping when composed with the inverse of any chart). Namely, 
in the coordinate system N, we have the following properties: 

(a) \gij{N{y)) - %| = o{\xy\) for all y G Reg{X); 

(b) I \yz\ - ||A^(i/)A^( 2 ;)|| 2 | = o{max{\xy\, \xz\})\yz\; 

(c) liniy^a; N(x)N{y)) = 0, where Ae stands for the standard Euclidean angle 

stands for the image through N of the direction of a geodesic from x to y; 

(d) \/.yxz —/.EN{y)N(x)N(z)\ = o{\yz\) when all the angles ofN(x)N{y)N(z) are bounded away 
from 0. (The same estimate holds if we consider a comparison triangle in and replace /.yxz 
by its Euclidean counterpart). 


Proof. In order to define N, we consider the bilinear form {gij)i,j defined in the chart F by the 
formula 


9io{z) = gij{x) + 


N „ 

k—1 


being g^j the components of the metric read in F and x = F{x), z = F{z). For z close to x, 
^ smooth Riemannian metric whose value at x coincides with that of {gij)ij and the 
same property is true for the first derivatives of gij. Therefore, we can find a smooth normal 
coordinate system for gij and we call N the coordinate system obtained by composing it with 
F. Note that in particular N is compatible with the DCq structure. By definition of a normal 
coordinate system, the metric components of g read in this normal coordinate system satisfy (we 
keep the same notation for simplicity) 


gij{x) = 5ij and 


^9ij 

dxk 


(x) = 0 for all i,j, k. 


Therefore, by construction, the same properties hold for g at N{x) when read in the chart N. The 
first item then immediately follows. Let us prove the second one. 

Set 7 a constant speed geodesic defined on [0,1] from y to z. Note that by density, it is sufficient 
to prove (b) for y, z G Reg{X). According to [TB], F{'j) belongs to C^((0,1)), therefore a := N{j) 
belongs to C^((0,1)) as well. More generally, Otsu and Shioya proved [161 Proposition 6.2] that 
the length L{9) of any curve 0 : [0,1] —>■ Reg{X) satisfies 

(50) L{0)=J Y^5e(s)(0'(s),0'(s))ds. 

Now, we infer from (a) the following estimate 


Ilk'(®)ll2-|a'(s)|,| 


||ct'(s)|| 2 + |cr'(s)|g 


< o(max{|a;y|, \xz})\yz\ 


where we use the fact that N is in particular biLipschitz (thus o(| • |) = o(|| • || 2 )). Integrating the 
inequality above leads to 


\\N{y)N{z )\\2 < A||.|| 2 (ct) < lyz] + o(max{|a:?/|, |a:z})|yz| 

where stands for the length of a with respect to Euclidean metric. The converse inequality 

is proved in the same way, up to replace the segment [N(y)N(z)] by an arbitrary close and parallel 
segment of same (Euclidean) length. Indeed, first notice that the Coarea inequality (see for instance 
[101 Theorem 2.10.25]) yields 


0 = jr^-^{Sing{X)r\R) > 


Jtf'^{Sing{X)r\p ^{z))dz, 
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being R a cylinder with axis [N{y)N{z)], and p the projection parallel to N{y)N{z). Therefore, 
by continuity, we can assume without loss of generality that fl Sing{X) = 0 so that 

the equality (PI) holds and the previous argument applies. 

We now prove (c). For any integer k > 0, set := 7 (^), being 7 a constant speed geodesic 
from X to 2 /, parameterized on [0,1] (so that toq = y, rui is a midpoint of x and y and so on). Note 
that for any vector u, 

lim lE{N{x)N{mk),u) = 

fc->-+oo 

Let us fix fc > 1. Elementary computations together with (b) lead to 

cos /-eN{ x)N{mk)Nirnk-i) + 1 < o(|a;TOfe_i|). 

Therefore, for y sufficiently close to x, we get 

^eN{ mk)N{x)N{mk-i) < \xmk-i\e{\xmk-i\) 

where limt_,.o £(t) = 0. By definition, |a;mfe_i| = ■ Moreover, up to reducing \xy\^ we can 

assume that |e(s)| < 1 on [0, |xj/|]. Consequently, we get 

00 -| 

fe =0 ^ 

and the proof of (c) is complete. 

It remains to prove (d). First, notice that since the angles of N{x)N{y)N{z) are assumed to 
be bounded away from 0, the error terms o(|fV(a;)A^(?/)|), o(| 7 V(a;)iV( 2 ;)|), o(|iV(z)A^(j/)|) are all the 
same and the property remains true if we use \xy\, \xz\, \yz\ instead since N is locally a biLipschitz 
map. Second, since these distances are assumed to be small, the law of cosines in the space form 
of curvature k implies that 

(51) \ cos Zyxz — cos Zoyxz\ < 0{\xy\'^) 

being /.oyxz the angle at a; of a comparison triangle in Now, using the Euclidean law of cosines 
for both N{x)N{y)N{z) and a Euclidean comparison triangle of xyz, and the estimate (2) (recall 
that o{\xy\) = o(|a;z|) = o{\yz\)) gives \ cos ZEN{y)N{x)N{z) — cosZoyxz] = o{\xy\‘^). Finally, 
combining this with (1511) yields 

I cosZyxz — cos ZEN{y)N{x)N{z)\ < 0{\xy\'^). 

Using again the fact that the angles of N{x)N{y)N{z) are bounded away from 0 (and tt as well) 
allows us to conclude. □ 

7.2.2. Proof of Provosition \7.5[ We will use the notation introduced in Proposition 17.81 By defi¬ 
nition, a DCq Riemannian manifold can be covered by countably many domains of (biLipschitz) 
charts, therefore it suffices to prove the result locally. Also, recall Theorem [TA] asserting that the 
metric g is differentiable .^^-almost everywhere (see Remark l7.7l) . Therefore, according to Propo¬ 
sition [7^ it suffices to prove the result at every point x G Reg{X) where g is differentiable and the 
first and second derivatives of / do exist, and we can use a normal coordinate system N around x 
to proceed. (Strictly speaking, to guarantee that it is indeed Hess“'^/ that appears in the Taylor 
expansion, we have to discard a .^^-negligible set of points where the second derivatives of / and 
the derivatives of the metric are not approximately continuous, and |IIess®/|(i?(a;, r) = o{r^) is 
not satisfied. For more details, we refer to the proof of Alexandrov’s theorem p242 in [3] or to 
the proof of Theorem 3.83 in [1].) For the rest of the proof, we fix such a point x and a normal 
coordinate system N around x. We have (writing / instead of / o A^”^) 

f{y) = fix) + dN^,)fiNix)N{y)) + iu^(,)/(iV(x)iV(y), iV(x)iV(y)) + o(|7V(x)iV(y)n. 
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Note that in a normal coordinate system around x, = Hessjv(a;)/- Therefore, in view of 

Proposition 17.81 we are done if we can prove 

(52) tN^f^],N{x)N{y) = oi\xy\). 

To this aim, we first claim that 

\^E -ZENiy)N{x)N{z)\ = |Z -ZEN{y)Nix)N{z)\ = o{\yz\) 

when all the angles of N{x)N{y)N{z) are bounded away from 0. 

The first equality follows from the definition of normal coordinate system.To prove the second 
one, take a point p such that N{p) is the plane N{y)N{x)N{z), N{x) is contained in the triangle 
N{p)N{y)N{z) and all the angles formed by these four points are bounded away from zero. Recall 
that for any triple x, y, z, one has Z. Zyxz by definition of Alexandrov space. Property 

(d) in Proposition 17.81 then gives 


^ \ X 1x5 lx 1x5 lx 

> ^eN {y)N {x)N {z) + ZeN{z)N{x)N{p) + ^eN {p)N {x)N {y) + o(|yz|) 
= 27r + o(|j/z|). 


On the other hand, Z txj'tx 27r (this is a consequence of the quadruple 

condition, see for instance 0 ) and the claim is proved. 

The proof of is by contradiction. Take a point yi such that \xyi \ = \xy\/2 and the direction 
between N{x)N{y) and N{x)N(yi). In particular, this gives 




> -lEN{y)N{x)N{y^) + lEf^^^,N{x)N{y) 

> o[\xy\) +l\%^^JyN{x)N(y). 


Thus if (IH^ were false, then we could construct a sequence yu ^ x with Z \^^^^^\N{x)N{yk) 
bounded away from 0 and this would contradict (c) in Proposition 17.81 


Appendix: The transformation law of the Christoffel symbols 


The goal of this part is to prove, in our nonsmooth setup, the following classical formula 
for the transformation of Christoffel symbols. 


(53) 


ab 


V dyo 




\m,t 


dXa dxb 


^ ^ d'^F, 


dye 


dXadXb 


17 —>■ R is a transition map and P, P 


used in the proof of Proposition 14.91 (see (l^rUl l. where F 
denote the Christoffel symbols in U and V respectively. 

We have not been able to locate a proof of (1551) in the literature which is only based on (155)) . 
So we provide such a proof for the sake of completeness. The tools we need to do so are basically 
those obtained in Corollarv l3.11l Let us recall that, by definition, 

dxj dxi dxi 


2n 


= E 


9 


For the sake of clarity, we shall use as much as possible matrix products to write down the formulas. 
For instance, the above formula can be rephrased in the following way 

dG\ ^kidgi. 

dxJ, 


2r^ = 


G 


-1 


dx^ 


ki 


kj 


E^ 


kl 

dxi ’ 


where G denotes the matrix ( 5 ^), i.e. the metric in the coordinate system of U. 
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We set Aj = G Denoting by G = {g)ij the metric in the coordinate system of V, we shall 

use the fact that (denoting by ‘M matrix transposition) 


G = *{dF){GoF){dF) 
and 

G-i = (G-ioF)‘(d;.(.)F-i). 


(54) 

Thus, 

(55) (^) =‘ (liridF)) (G o F) (dF) +‘ (dF) {^{G o F)) (dF) +‘ (dF) (G o F) {^{dF)). 




'dxi 


'dxi 


Accordingly, we decompose Aj into three terms: 


A,- - An + Bn + Gn, 


where 


(56) 


'A = {dFi.)F-^){G-^oFY{dp^.^F-^y{£-idF)){GoF){dF) 

B, = {dp^.^F-Y{G-^oF){£-{GoF)){dF) 


Cj = {dF^.)F-Y{^idF)). 

We first consider Gj and we compute {Cj)ki- 


(c,H. = y:^oF 


dys 


dxndxi 


Using the symmetry of the second distributional derivatives, we get 

Now, we treat the term Bj. We shall use that for h G BVq (see m in Corollary 13.111) . the 
chain rule formula holds in the following form 


(57) 

Thus, 


!-('-) 4 


= E ^ [G(.)r-') (G-‘ o F) (F-(^)){dF) 


and we get by expanding this expression 


m.s.t 


dxj dys 


(G-‘or)(f.(|£)) 


Similarly, 


{Bi)kj — 




dFt dF^ 

dxi dys 


-1 


oF 


(g-^ o f ) 


By adding the two formulas above, we get 


{Bj)ki + {Bi)kj — 


dFm dFt dF-^ 

— dxj dxi dys 

m,s,t d 


g-^of) if* 

dym ) 


- ki 

dFt 
dxi' 

dFrr. 


dxj ■ 


oFF* 


G- 


+ G 


_^dG 
dyt 


where we use the fact that for a € BVq, g € QA4, F*{afi) = {ao F) F*{fj,) (see Corollary 13.11|) . 
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Let US notice that 


dG 




dV'. 


dy, 


= F‘ 


(2r^*) +E6 


dye 


gse oFF* 


dye 


Therefore, if we combine all the equalities above, we get 


2(r'=.),.=2^^oF 


O dF^ dFt dF^ \ u u* /"f* ^ 


dys dx-jdxi ^ dxj dxi dy 

J nm <i i 


y dF^ dFt dF^^ 

m^t,e 


oFg‘^^ oFF* 


dyrnt 

dye 




dxi. 


The proof is then complete if we prove that 


(58) E 


dFra dFt dF^ 

m,s,t,e ^y^ 


-1 


dgv 


„ F J.% ( ^ ) + + {/!.)., -T.S‘'£^ = 0^ 


dxu 


To this aim, we write according to 


iAj)k^= E 


dys dyt dxjdxu 


dxi ' 


By writing {Ai)kj the same way (exchanging the role of v and w) and then adding the two 
expressions, we get 


(59) + {A,)k^ = E ^^°Fr*oFg,^oF^^oFj^(^^). 

dvs dyt dxu V dxi dxi) 


dFz^ ^ d fdF^dF, 


dy. 


dyt 


To conclude, it remains to rewrite 'Ylii9^^F^- Using (IH^ . we obtain 


dgij _ d ,dFyj dFy^^ dF^, dFy d , 


dx. 


dxi dxj dxu 


Then, combining this together with (15^ : 


^ku 


dF, 


= V -T^oFg^^ oF 
^ dy.^ 


S,t 


dF- 

dyt 


oF 


yield. 


, 60 ) i: :^.FA^oF^(.^'j^%,oFroF 


dF, 


dF- 


d , dFy, dFy 


dxu ^ , dys 

u,w,v,s.,t 


dyt 


dxu ^ dxi dx-j 


(D) 


^oF^oF^^^igy,yoF)roF. 
, dys dyt dxi dxj dxy 


(B) 
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According to (15^ . {Ai)kj + {Aj)ki = {D). To complete the proof, we apply (1571) to g^v in [E). 
This gives 

{E)= Y. 


= E 


dvs 


-^°F 

avs 


dyt dxu dxi dxj ^ dye 


dFw dFy ^ dywv \-^s9 
dxi dx 


t 


F*{^)roF. 

dye 


Therefore (E) coincides with the first term in (15^ and the proof is complete. 
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